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Abstract 

We provide multiple Schramm-Loewner evolutions (SLE) to describe the scaling limit 
of multiple interfaces in critical lattice models possessing Lie algebra symmetries. The 
critical behavior of the models is described by Wess-Zumino-Witten (WZW) models. 
' r-| \ Introducing a multiple Brownian motion on a Lie group as well as that on the real line, 

O j we construct the multiple SLE with additional Lie algebra symmetries. The connection 

between the resultant SLE and the WZW model can be understood via SLE martingales 
satisfied by the correlation functions in the WZW model. Due to interactions among 
SLE traces, these Brownian motions have drift terms which are determined by partition 
functions for the corresponding WZW model. As a concrete example, we apply the 
formula to the su(2)fc-WZW model. Utilizing the fusion rules in the model, we conjecture 
that there exists a one-to-one correspondence between the partition functions and the 
topologically inequivalent configurations of the SLE traces. Furthermore, solving the 
Knizhnik-Zamolodchikov equation, we exactly compute the probabilities of occurrence 
for certain configurations (i.e. crossing probabilities) of traces for the triple SLE. 

1 Introduction 

Geometric aspects of critical phenomena are characterized by random fractals such as con- 
formally invariant fluctuations of local order parameters. They have been extensively studied 
from various different points of view, especially in two dimensions (2D) where the conformal 
invariance imposes strong constraints on the structure of critical phenomena. Among them, 
the Schramm-Loewner evolutions (SLEs) pQ, which directly describe geometric aspects of 2D 
critical phenomena through simple ID Brownian motions, have brought a renewed interest 
in the theory of random fractals (see [21 [31 HI El El El El [9] for reviews). 

The SLE is a stochastic process defined in the upper half plane H. Its evolution is 
described by the ordinary differential equation 

dg t (z) = —^- — , go(z) = zei, (1.1) 
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Figure 1: Schematic view of a SLE trace (74 denotes the tip of the trace) (a) and its hull Kt 
(b). 

where xt = v^£t * s a Brownian motion on M, starting at the origin (i.e. xq = 0), and its 
expectation value and variance are given by E[da?t] = and ~E[dxtdxt] = ndt, respectively. 
Here k > is a diffusion coefficient which essentially characterizes the SLE process. The 
SLE (jl.lj) has a solution up to the explosion time t z , i.e. the first time when gt(z) hits the 
singularity x t . Let K t = {z £ H|r z < t} be the hull at time t (see Fig. Q] for a schematic 
view). Then Kt (t > 0) is an increasing family of hulls: K s C Kt for s < t. Moreover gt(z) 
with gt(z) = z + 0(l) at z — > 00 (hydrodynamic normalization) is the unique conformal map 
uniformizing the complement of the hull Kt in the upper half plane H: gt(z) : H \ Kt — > H 
(see Fig. [2]). The image xt by g^j~ (xt) defines the tip 74 of the growing random curve. More 
precisely, it is expressed as 74 = lim^+o g^~ 1 (xt + ie). 

The connection between the SLE (jl.ip and conformal field theory (CFT) is well under- 
stood [IHlIIIlIElIiaElIiaillKIZlElE]. Specifically, it can be accomplished by noticing 
that CFT correlation functions 

Mt = MocWfa)) (L2) 

are SLE martingales: dMt/dt = 0, where ipilt) and ip(oo) are boundary condition changing 
(bcc) operators with conformal weights h, which are inserted at the points z = jt and z = 00, 
respectively (see |T7] or next section for details). Thus one can find the SLE corresponds 
to the minimal conformal field theory M(p,p') (p, p' are coprime integers satisfying the 
condition p > p' > 2) where the central charge and the conformal weights of the primary 
fields are, respectively, given by [181 EH] 

n c (P - v'f , far ~ P's) 2 - (p - v'f n „v 

c = l-6 — , h rs = — . (1.3) 

pp' 4pp' 

Then the diffusion constant k and the conformal weight h of the boundary field if) are, 
respectively, expressed as 

c _ (3/c-8)(6-/c) ^ _ 6 - k _ (h 2 ,i for K = Ap'/p< 4 
2k 2k for k = 4p/p' > 4 

Namely the boundary field ip is degenerate at level two, i.e. possesses a null field at level two. 



2 




Figure 2: Uniformizing map g^. EI \ Kf — > M and its inverse. 

An extension to the SLE connecting with conformal field theories with Lie algebra sym- 
metries (i.e. WZW models [2Q|. |2T| I22j ) has been achieved by adding the extra Brownian 
motion exp(^ a t a d8f (z)) on a (semisimple) group manifold G associated with a Lie algebra 
g, where t a, s (a = 1, . . . , dimg) stand for any representation of the Lie algebra generators 
[23} [23] (see [25] for a very different approach). The evolution of this additional stochastic 
process Of (z) is defined as 

d9?{z) = E[fM a] = 0) E[d-d?d$] = 5 ab dt {r t G R). (1.5) 

Z Xt 

The combination of (jl.ip and (ll.5p defines a fractal curve living on the Lie group G manifold. 
The SLE martingale dAit/dt = 0, which should be satisfied by the CFT correlation function 
(|1.2p (note that the bcc operators and the operator O take their values on G) , determines the 
relation of the SLE defined by (jl.ip and (|1.5p with the corresponding WZW model. Some 
extensions to the SLE with other additional symmetries have also been done in |26j ¥27\ |2"81 

1221 EDI EQ. 

In this paper we generalize the SLE for the system containing multiple random interfaces 
which possess additional Lie algebra symmetries, according to the theory developed in [17] , 
Assuming that each SLE interface grows under an independent martingale in the infinitesimal 
time interval, we extend (|l.ip together with (jl.5p to the case for the system with multiple 
interfaces. The evolution (cf. (jl.ip for the single case) characterizing the geometric aspect of 
the interfaces is described by a multiple Brownian motion on the real line, while the evolution 
(cf. (jl.5p for the single case) expressing the algebraic aspect is described by a multiple 
Brownian motion on the Lie group G. Both the Brownian motions, however, have drift terms 
describing the interaction among the SLE traces. Taking into account the SLE martingale, 
one finds that these drift terms are determined by the partition function in the corresponding 
WZW model. Moreover this partition function characterizes the configuration of the SLE 
traces. As an example, we apply our formula to the SLE for the su(2)fc-WZW models. 
Utilizing the fusion rules for su(2)fc, we conjecture that there is a one-to-one correspondence 
between the partition functions and the topologically inequivalent configurations of the SLE 
traces. Further, we exactly compute the probabilities of occurrence for certain configurations 
of traces (that correspond to crossing probabilities) for the triple SLE, which can be obtained 
by solving the Knizhnik-Zamolodchikov (KZ) equation |22j . 

The paper is organized as follows. In the subsequent section, we describe some basic 
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notions required in this paper. In section 3, the multiple SLE for conformal field theory with 
Lie algebra symmetries is formulated. The drift terms of the driving Brownian motions are 
explicitly determined in section 4. A concrete application of the formula to the SLE for the 
su(2)fc-WZW model is given in section 5. 



2 Preliminaries 

In this section, we describe several theoretical foundations required in subsequent sections. 
In the former part of this section, we introduce SLE martingales from the point of view of 
statistical mechanics (see [7] for details). SLE martingales are given by CFT correlation 
functions involving bcc operators, and hence they become a key to decipher the relation 
between SLE and CFT. In the latter part, general properties of correlation functions are 
explained in the case of the WZW model. 



2.1 SLE martingales and conformal correlators 

To formulate the SLEs correctly describing the behavior of 2D interfaces in critical systems, 
one must construct SLE martingales in terms of the corresponding statistical systems defined 
on H. Let (O) be the thermal average of an observable O defined in H, and (0)|{ 7t } be 
the thermal average under a given shape of configuration of (multiple) interfaces, where {7^} 
denotes the shape of configuration with its occurrence probability P[{7t}]. Then the thermal 
average (O) must be given by 

(O) = E[(0)\ H} ] = J2 P li^}](0)\ M , (2.1) 

M 

where the average E[- • • ] should be taken over all the possible configurations. The conditional 
expectation value (C)|r 7 j is thus time independent (i.e. conserved in mean), and therefore it 
is a martingale. Here and in what follows, we denote it by Mf 

At the critical point where the conformal invariance is expected in the system, the above 
observable can be described in terms of the CFT correlation functions. For the situation 
where the number of the interfaces under consideration is m, it reads 

where Wj's denote the positions of the tips of the interfaces. Note that the CFT correlation 
functions are defined on the domain M removing the hull Kt, i.e. M\ Kt (see Fig. [3]). The 
operators ^,-(u;,-)'s inserted at the positions of the tips denote primary bcc operators with 
conformal weights h: under a local conformal map z — > z' = w(z), a primary field ip{z) 
transforms as ip(z) — > ip'(z'): 



dw(z) 



-h 



vco = y-ir) (2-3) 

The denominator in (|2.2[) stands for the CFT partition function with a specific boundary 
condition fixed by the bcc operators ipj(wj). 
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(C'0l(Wl) ■ ■ ■ lp m (w m )lp m+1 {oo)) { gi 04'i(x U ) ■ ■ ■ i! m {x m t)lpm+l{oo)) 

(^l(wi) • ■ ■ 1p m (w m )ll} m+1 (oo)} (Mxit) 1pm{Xmt)i'm+l(oo)) 



Figure 3: CFT correlation functions on EI \ IK (a). It can be transformed to the one defined 
on EI by the uniformizing map gt(z) (b). 



Applying the conformal uniformizing map gt(z) (written with the same symbol as that 
used for the single SLE (jl.ip ). we obtain 

t (4>l{xu) ■ ■ ■ i>m{x m t)i>m+\{oo)) 

Here Xjt = gt(wj), and 9t O is the image of O by the map gt- Note that the Jacobians coming 
from the conformal map on ipj have been canceled in the numerators and the denominator!^]. 
Now SLE martingales Ait are expressed as the CFT correlation functions on EI, where the 
bcc operators ^j(xj) are inserted at the points Xj G K (see Fig. [3]). 

To proceed further, let us consider the case when the operator O is a product of an arbi- 
trary number of primary fields O = TYj^&jizjjZj) a t positions (zj,Zj) and with conformal 
weights (hj,hj). By construction, the uniformizing map gt{z) can be analytically extended 
to the lower half plane: gt{z) = g~t(z). Then, the doubling trick can be applied to the CFT 
correlation functions. The result reads 

^ _tt f^Vjt\ hj <nj=i Myjt)M x u) ■ ■ ■ Vv»(xmt)Vwi(oo)) 

4 J=1 V^J {lpl{xit) ••• i>m(x m t)ll>m+l{oo)) 

where we denote that y jt = g t (zj); 4>j +n = 4>j] z j+N = zj = z*; y j+ m = Vjt = Vj t , h j+N = hj 
for 1 < j < n. Here <pj(yjt) {4>j{Vjt)) stands for the holomorphic (antiholomorphic) part of 
the field (f>j(yjt,y~jt)- 

In this paper, we analyze the SLE martingales (|2.5[) for the system that possesses addi- 
tional Lie algebra symmetries. Namely we construct the multiple SLE for the WZW models 
which are one of the most fundamental CFTs with extra Lie algebra symmetries. In this case, 
the primary fields constructing the SLE martingale (|2.5p possess internal degrees of freedom, 
such as "spin". 



1 The identity (z[) ■ ■ ■ &' n (z' n )} — {$i(z' 1 ) ■ ■ ■ & n (z' n )) holds for a global conformal map z — >■ z' — f(z), 
where the field transforms as Qj(zj) — > &j(z'j) [19] . 
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2.2 WZW models and correlation functions 

Let us introduce several properties for correlation functions of WZW primary fields to analyze 
the SLE martingales (|2.5p . The WZW model is a CFT described by a field g(z,z) taking 
values in a group manifold G associated with a Lie algebra g [X9^ |20~1 \21\ [22] . The model is 
invariant under 

g(z,z) -+n(z)g(z,z)Q- 1 (z), (2.6) 

where and Q(z) denote arbitrary matrices valued in G. This invariance gives rise to Noether 
currents J(z) = —kd z gg~ l and J(z) = kg~ l dzg, which can be written as 

dimg dimg 

= E ja (^ a < J >~) = E Ja (^ ( 2 - 7 ) 

a=l a=l 

where i a, s stand for any matrix representation of the generator of q, with commutation 
relations [£ a ,t 6 ] = ^ c i/^f. The parameter & is a positive integer referred to as the level. 
Hereafter we only consider the holomorphic components, as if there were no boundary (cf. 
(|2.5p ) [32j. Let {X) be a correlator of G-valued fields. Then the infinitesimal transformation 
£l(z) = 1 + uj{z) leads to the Ward identity 

„ dim g 

S U (X) = -— fdzY, u; a (z)(J a (z)X). (2.8) 

m J a=l 

On the other hand, conformal aspects of the WZW model are described by the stress 
energy tensor T(z): 

dim g 
a=l 

where h v denotes the dual Coxeter number of g. Note that : : means the "normal ordering" 
defined as 

1 f dw 



: A(z)B(z) : = <f> A(w)B(z). (2.10) 

2-rn J z w — z 

The infinitesimal conformal transformation z — > z' = z + e{z) leads to the Ward identity 

S e {X) = -^-jdze{z){T{z)X). (2.11) 

This geometric part of the Ward identity (|2.1ip together with the algebraic part (|2.8p are a 
key ingredient in analyzing the SLE martingale. 

The current J a {z) and the stress energy tensor T(z) can be expanded in terms of the 
modes .1% and L n , respectively. Namely 

J a (z) = E J>- n -\ T(z) = E L n z~ n ~\ (2.12) 
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where J% and L n denote, respectively, the generator^ of the affine Lie algebra g with level 
k, and the generators of the Virasoro algebra: 

c 

[L n , L m ] = (n - m)L n+m + j^(n 3 - n)5 n+mfi . (2.13) 



The central charge c of the Virasoro algebra is then given by 

k dimg 



(2.14) 



By construction (|2.9|) . the Virasoro generators are not independent from the affine generators: 

dimg 

Ln = 0(h _l hV\ 1 Jm^n-m : > ( 2 - 15 ) 

z{k + n ) a=i mgz 

where the normal ordering : : means that the operator with larger index n is placed at the 
rightmost position. 

The primary fields in the WZW model are defined as the fields transforming covariantly 
with respect to the G{z) transformation ()2.6p : 5^9 = wg. Together with the conformal 
covariance (see (|2.3p ). these properties can be expressed in terms of the operator product 
expansions (OPE) via the Ward identities (|2.1ip and ([2.80 : 

T(z)ipx{w) = 7 v7 H ^ r eg-, J (z)^ x (w) = — - hreg., (2.16) 

(z — w) A z — w z — w 

where the field ijj\(z) takes values in the representation specified by the highest weight A, 
and tfi is the generator t a in that representation. Furthermore utilizing the field-state corre- 
spondence, i.e. |V>a) : = hniz-K) ^a(^)|0), we can translate these properties into 

L \^x) = /ia|^a), L n \^x) = (n > 0), 

■Wa> = -Wa>, «a> = (n > 0). (2.17) 
Insertion of the relation (|2.15p into the l.h.s of the first equation in the above yields 

1 dimg dimg . , 

v ' o=l a=l v ' 

In the last equality we used the explicit form of the eigenvalue of the quadratic Casimir. The 
quantity p denotes the Weyl vector, i.e. the sum of the fundamental weights p = Yli=i A*> 
where r is the rank of g. Comparing the r.h.s in the first equation in (|2.17p with the above, 
we find 



2 Here we use the orthonormal basis in terms of the Killing form: K(J%, J b n ) — 5 n + m fiS a ' b . In that case, 
the structure constants can be written as f ab c — f a b c , where f a bc is antisymmetric in all three indices. 
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Since the Virasoro generators are expressed in terms of affine generators as (|2.15p . the 
other arbitrary states are of the form «7*j J_ t ■ ■ ■ with l\,l2 ••■ positive integers. Let 
(L_iip\)(z) and (J^_ l i/)\)(z) be the descendent fields corresponding to the states L_i\ip\) and 
J-il^x), respectively. Combining the OPE (|2.16p and the mode expansions (|2.12p . one finds 
that the correlation functions ((L_iip\)(z)X) and {(J® l ip\)(z)X), where X = n^=i V'AtO^)) 
satisfy the following equations: 



{(L^ x )(z)X) = c^ x {z)X), C-i = E 



3=1 



(I - l)h Xi d; 



( Zj - z y ( Zj - z y-i 



A i 



((j a _M(z)x)=j» l (Mz)x), ^ = Ett^v- ( 2 - 2 o) 



i \ Z 3 ~ Z > 
3=1 J 



Note that the global G-invariance requires 5 L0 (tp\(z)X) = 0. Using the Ward identity (|2.8|) 
with constant u and the OPE (|2.16p . we obtain a relation satisfied by the correlation function: 



dim t 



E K + E^ (M*W=0- (2-21) 

0=1 \ J=l / 

This constraint together with the global conformal invariance (or equivalently SL(2, C)- 
invariance) fix the structure of the two- and three-points correlation functions. 

To close this section, let us derive a crucial equation called Knizhnik-Zamolodchikov 
(KZ) equation |22] satisfied by the correlation functions of the WZW primary fields. The 
constraint stems from the fact that the Virasoro generators are not independent from the 
affine generators as in (|2. 15|) . For n = —1, we have L_i = Y^ a J-i^o/i^ + ^ V )- Then a null 
state \x) = is given by 

( dim \ 
L-l + TT ^ 7 E *S J -i) M = °> (2.22) 

where we used the property in (|2.17p . Using the field-state correspondence and inserting the 
property ()2.20p into the correlation functions (x(z)X), one obtains the KZ equation: 

dimg m f a f a \ 



o=l j=l 



Here we used the translation invariance (d z +Y2j 9 Z j)(x( z )X) = 0, which can be easily verified 
from (l2Tlj) and the OPE (12TTBJ) by setting e(z) = e. 



3 Multiple SLEs for WZW models 

Now we generalize the SLE (jl.ip and (|1.5p for the system containing multiple random inter- 
faces with additional Lie algebra symmetries. In the infinitesimal time interval, we expect 
that each SLE interface grows under an independent martingale. Let us discuss the case where 
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the number of the interfaces is m. Then the uniformizing map gt(z) with the hydrodynamic 
normalization may be of the form [T7j (see also [33 } 134 4 [35] for other approaches). 

m ^ 

d 9t(z) = V . r — , g {z) = z, (3.1) 

^gtiZj-Xat 

where dq^s mean infinitesimal time intervals satisfying the condition YJa=l ^la = dt. The 
random processes x a (1 < a < m), which play a role as driving forces for the growth of 
interfaces, should be written as the Ito stochastic differential equations: 

dx at = y/ndZat + dF at , (3.2) 

where £ a t is an R m -valued Brownian motion whose expectation value and variance are, re- 
spectively, given by 

E[d&rt] = 0, E[dUd&t) = $ a pdq a . (3.3) 

Namely dq a prescribes the growth rate of each interface. The quantity dF a t denotes a drift 
term proportional to dq a , which comes from interactions among interfaces, and will be de- 
termined later by the SLE martingale. 

To extend (|3.ip to evolutions with Lie algebra symmetries, we define a stochastic process 
ex PEa d9f(z)t a ] (a = l,...,dimg) living on a Lie group manifold G(z), where 6f(z) is 
written as 

m rl a 

dB ?W = E7??-' TO = (3.4) 

z x a t 

a=l 

with 

dp at = ^dO a at + dG at . (3.5) 
Here is an M. m dim -valued Brownian motion with 

E[d$ a at ] = 0, E[dKtd$U = S ab 5 a(3 dq a , (3.6) 

and dG at stands for a drift term proportional to dq a , which will be also determined later. 
For z = xpt, we must define 

m A a 

de a t {Ht) = Y. • ( 3 - 7 ) 

i X/3t- X a t 
a=l H 

The growth of interfaces with affine Lie algebra symmetries is described by both the 
geometric (|3.1|) and the algebraic (|3.4|) components. 

4 Drift terms and SLE martingales 

The remaining problem is to determine the drift terms appearing in the driving forces (|3.2p 
and in the Brownian motion (|3.5p . To achieve it, we must evaluate the variation of the SLE 
martingale A4t involving WZW primary fields. For the WZW models, the SLE martingale is 
written as 

Mt = f\ ( 9yjt ) ^ ^j=i <Kj (jftOjfri (xu) ■ ■ ■ ^x m (x mt ^ Xm+1 (oo)) 
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where yjt = gt(zj) and x a t = gt(w a ) (see section 2.1 for details). Note that the primary fields 
(fipj and i/j\ a constructing the correlation function take values in the representation specified 
by the highest weights fJLj and X a , respectively. 

The drift terms are determined by the condition which makes Mt to be a martingale. 
To simplify the notations we sometimes omit the index t, A and [i (e.g. x a t = x a , 9f = 
® a i ^Aa(^a) = ^aixa), h\, = hj, etc.), if there is no ambiguity. Let Z t , zf and jf be 
the denominator, numerator and Jacobian factor of (|4.ip . respectively. First we explicitly 
compute djf . A simple manipulation leads to 



M'sr*(^)-*(i)^(s 



2h 



dy 

dzj (y-x a )' 



r, (4.2) 



where in the last equality, we applied the geometric component of the SLE (|3.ip . Thus we 
obtain the variation of the Jacobian factor: 



m 2n 



dJ? 



2flA 



(4.3) 



a=l j=l v ^ 

Next we shall calculate d(jfzf ). The ltd derivative of the bulk fields <f)(y) is given by 

dimg dimg dimg 

2 



My) 



d(f>(y) dy 



E fr dq « + E ^dMv) + J E * adr M E ^(vMv) 



dy ^dq, 



a=l 



a=l 



6=1 



E 

o=l 



^ y - x Q 2 ^ (y - :r a )< 



y - x c 



a=l 



0(y). 



(4.4) 



Here we applied (|3.ip to the first term, and (j3.4[) to the second and third terms. For the third 
term we also used the property (|3.6j) . Similarly, for the boundary fields ip a (x a ), we obtain 



dim g 



dim g 



dim g 



dq a ^d 2 Xa + dx a d Xa + £ t a d^(x a ) + 1 jr » (*«) £ t fe d^( 



a=l 

dim g ra 



a=l 



6=1 



0a (x c 



(4.5) 



The above two relations together with (|4.3p give the form of d{jf Zf). Explicitly it reads 

dim g \ 



Mm 
jf 



a=l 



j^a,-l ~ 2 £a,-2 + 2 E ^a-l^a,-l \ Z t 
dim g \ 



+ 2C a - 2 dq a + £ a ,-idx a + Ja,-idp a a Zf, 



(4.6) 



a=l 



a=l 



10 



where we define 



_ , ( Z _ 1)hj _ d yj \ i)h p a 



ja V 7 ^t + V ? ^-r, (4.7) 

and C (resp. J) is given by subtracting the first sum depending on yj from the r.h.s of the first 
(resp. second) equation in (|4.7p . Note that the operator C (resp. J) characteristically appears 
in the correlation functions between a descendent field (L_iip\ a )(x a ) (resp. (J_it/j\ a )(x a )) 
and some composite primary field (see (|2.20p ). To derive (|4.6p . we have also applied the 
identity d Xa Zf = C a -iZf, which comes from the translation invariance of the correlation 
functions (see the explanation below <|2.23[) ) . 

In completely the same manner, dZt can also be evaluated: 

m / dimg \ 

dZ t =Y,dq a [ ~ 2£ Q ,-2 + ^a,-Ja,-l Zt 

0=1 V a=l J 

m / dimg \ 

+ 2£*,-2dq a + C a ,-idx a + Ja,-idp a a Z t , (4.8) 

a=l \ a=l J 

By construction, we must set h\ a = h (1 < a < m) for the conformal weights of the 
bcc operators. Furthermore, due to the constraint between the conformal weight h\ and the 
representation A (|2.19p . all the bcc operators ip\ a (x a ) except for tp\ m+1 (oo) must take values 
in the representation specified by A, or its conjugate A* (note that h\ = h\» holds). The 
conformal weight h\ m+1 for the bcc operator inserted at 00 is determined by fusion rules. 
Thanks to this restriction, we can determine, in principle, the relation of the parameters k, 
t and the weight h\ of the bcc operators by considering the case that there is only a single 
interface (m = 1) as in |23|, [24"] . We see this in the next subsection. 



4.1 Single case (m = 1) 

For the single SLE (set a = 1 and m = 1), one sees dZt = due to the translation invariance 
of the two-point correlation functions (it can also be obtained directly from (|4,8p by setting 
m = 1). Therefore the SLE martingale gives a constraint to the correlation function Zt, i.e. 
~E[d{jf zf)] = 0. From the relation (|4.6p and the definitions (|3.2p and (|3.5p . one easily see 
that this restriction leads to dF\ = 0, dGf = and 

(dimg \ 
1 Ll_ x - 2£ a ,_ 2 + l£ J^-iJZ-i J Z* = (a = l). (4.9) 

Namely any drift terms do not show up in the driving forces (|3.2p and (|3.5p . as expected. 
Moreover the constraint (j4.9|) indicates that the bcc operators must have null states at level 
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2. Namely Mt is a martingale if and only if the bcc operators have the following null states 
at level 2: 

dim g 



0=|X> ■- 



(dimg \ 



(4.10) 



Here L n and J% are, respectively, the Virasoro and affine generators (|2.13[) . This is equivalent 
to the conditions J\\x) = and J||x) = 0, which respectively lead to 



(rife - r/i v - 2)J b _ x + KJ b L_! + ir f ab c J$J-! 



0, (« + r/i V -4)J 6 | 



0. 



J4.ll) 

This necessary and sufficient condition is rather involved. For some simple cases (e.g. iu(2)^), 
however, we can directly solve the above equations by acting the generator jf . (Note that 
more elegant procedure utilizing the KZ-equation has been developed in [24 J). 

For later convenience, here we explicitly write down the results for the iu(2)fc-WZW 
model. The central charge c (|2.14|) and the conformal weight h\ 1 (I2.19|) of the bcc operator 
r tjj\ 1 (xi) in the spin-A/2 representation are, respectively, written as 



3k 
k + 2' 



A(A + 2) 
4(fc + 2) ' 



(4.12) 



where A denotes the fundamental weight, and we used h v = 2. From the direct evaluation of 
(|4.1ip for su(2)/% case, one finds that the conditions in (|4.1ip are valid only for the case that 
the bcc operator carries spin-1/2 (A = 1) [24] . Then [23} I24j 



4(fc + 2) 



(jfc > 2). 



(4.13) 



fe+3 ' k+3 

For k = 1, k and r can not be specified, and only the relation k + 2r = 4 is imposed. This 
case, however, corresponds to a c = 1 CFT (cf. (|1.4|) ), and therefore we shall set k = 4. 



4.2 Multiple case (m > 1) 

Now we identified the bcc operators ip\.(zj). Namely they have null states at level 2 and 
must satisfy the condition (|4.10p . Utilizing the field-state correspondence, one finds that the 
first sums in (|4.6|) and (|4.8p vanish. 

Thus the drift term of the Ito derivative of the CFT correlation function 



dM t = d 



jjzl 
z t 



d{J?Z?) J?zfdZ t d{J?Z?)dZ t , JfZXdZ, 



<t> ryfyl 



zl 



zl 



+ 



z? 



(4.14) 



is explicitly given by 

m 

V[dM t ] =4 E dF « ~ Kdq a d Xa log Z t - 2 J2 ■ 

13=1 
/3^a 



a=l 



dq p 



X0 



0, 




+ 



J 



m dimg 

EE 

a=l a=l 



dG° 



Tdq a ~ a 
Z t 



Jn, -\Zt 



ja y<P 



z*j- 



a 

a,-l 



z t 



Zi 



(4.15) 
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where we substituted the relations (|3.2p and (|3.5p . By recalling that Mt is the SLE mar- 
tingale, the above quantity must be zero. Thus one finds the drift terms dF a and dG a are 
described by the partition function Zt: 



dF a = ndq a d Xa log Z t + 2 ^ - 



A/3 1 



T 



dq a . 



(4.16) 



4.3 Main claim 

To summarize, we have constructed the multiple SLEs for g^-WZW models: 

2dq a 



d 9t(z) = 

a=l 
rn 



{ 9t{z) - x a t 
dplt 



a=l 



J at 
Z Xat 



dx a t = V^d^at + dF at 
dp a at = ^dr at + dG a at (1 < a < dim(fl)). 



(4.17) 



The drift terms dF a t and dG a at are given by 



dF^ = Kdq a d Xat log Z ( + 2^ 



dGat 



m j.a 7 

Z * ^1 X ^ ~ Xat 
P^a 



(4.18) 



where £ a t (resp. is the R m (resp. R mdim9 )-valued Brownian motion whose expectation 
value and variance are given by 



E[dCot] = 0, E[d£ at d€i3t\ = 5 a/3 dq a , 
E[dd%t] = 0, E[d# a at d#%] = 5 ab 5 a pdq a . 

Zt is a partition function involving the bcc operators: 

Zt = (tpXi (xu) ■ ■ ■ ip\ m (x mt )ip\ m+1 (oo)) , 



(4.19) 



(4.20) 



where all the bcc operators ip\ a (x a ) except for ip\ m+1 (oo) must take values in representations 
specified by A, or its conjugate A*. The conformal weight h\ m+1 for the bcc operator inserted 
at oo is determined by fusion rules (see next section for example). The structure of the 
partition function is described by both the global G-invariance (cf. (|2.2ip ). and the KZ- 
equation (cf. (12331) ) 



dimg m+1 
a=l a=l 



( 



7 



dimg m -^a 



\ 



i ~ — ° ■■- 

EE 



k + h v 



a=1 p ^a-Zp 

P+o. 



0. 



(4.21) 
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5 Multiple SLEs for iu(2) fc -WZW models 



As a concrete application of our formula, let us consider multiple SLEs for the su(2)fc-WZW 
model where the bcc operators carry spin-1/2, and discuss topologies for the SLE traces. 

For the su(2)fc case, the fusion rules of the primary fields are similar to those in the 
minimal CFTs. Therefore it is natural to extend the argument [T7] describing the geometric 
configurations of the SLE traces in minimal CFTs to the iu(2)fc-WZW model. Thus we make 
the following conjecture. 

Conjecture 5.1 There exists a one-to-one correspondence between topologically inequivalent 
configurations o/su(2)fc multiple SLE traces and the independent solutions of the KZ-equation 
satisfied by the partition functions. 

We describe this conjecture more specifically. Let x a (1 < a < m) be the positions where 
the m SLE traces start to grow, and be ordered as x\ < ■ ■ ■ < x m < oo. Consider the case 
that the traces eventually form m — n disjoint curves in H so that each point x a is an end 
point of exactly one curve and oo is an end point of exactly m — 2n > curves. Namely n 
disjoint curves form arches and other m — In curves converge toward the point at oo (see 
Fig. [3] for m = 4 and n = 2). Then the number of topologically inequivalent configurations 
are given by 




This is nothing but a Kostka number appearing as the coefficient of the irreducible decompo- 
sition for the m tensor product of the su(2) fundamental representation L® m into £( m _2n)A 
(Laa denotes the integral representation with the highest weight AA): 

\m/2\ 

•^A = c m,vX'm—2ni (5-2) 

n=0 

where |_^J denotes the integer part of x. 

Now we mention that the relation between the geometric configurations and the CFT 
partition functions Z (|4.20[) (hereafter we omit the index t to simplify the notation). To this 
end, we consider the fusions of the su(2)& primary fields: 

min(A+/j,2fc— A— fj.) 

ip\A ® VVA = (J) VVA, (5.3) 

u=\X— fj,\ 
A+/j+!/=0 mod 2 

where ip\\ stands for the su(2)/<. bcc primary field taking values in the integral representation 
Laa- For sufficiently large k, the above rules reduce to the decomposition (j5.2|) . Applying the 
fusion procedures (|5.3p to the bcc operators ij)x a 0- — a — m ) recursively, one can reduce Z to 
a two point function involving the fusion operator and the bcc operator at oo. Thus, for the 
non- vanishing partition functions, the conformal weight h\ m+1 of the bcc operator ip\ m+1 (oo) 
must be equivalent to that of the fusion operator. Fig. [5] shows the possible conformal weight 
h\ m+1 up to m = 4. Due to the fusion procedures (|5.3p . for sufficiently large k, one finds that 
the number of the paths from the left to hr m _ 2r AA corresponds to the number of topologically 
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distinct configurations for the SLE traces, i.e. c m ^ n . For generic k, the number of paths is, in 
general, constrained by the fusion rules, which affects the structure of the partition functions, 
and therefore the realization of the geometric configurations. 

For instance, the configuration where all the curves eventually converge toward oo (i.e. 
no arches) corresponds to the path to the weight h\ m+1 = h m \. By a standard argument for 
the CFT correlation functions, one obtains the corresponding partition function Z: 

m 

Z = \\{xj - x fc )2p+2) for m<k. (5.4) 

j<k 

One can easily check that this partition function satisfies the KZ-equation ()4.2ip . The fac- 
torized correlation function with the same exponents \/{2{k + 2)) does not exist for m > k, 
indicating that the no-arch configuration is allowed only when m < k. 



5.1 Double SLEs 

According to the argument in [17j . we confirm our conjecture by considering the simplest case 
(m = 2), where only two SLE traces exist in the upper half plane HI. To analyze specifically, 
let us write down the partition function Z. Up to a constant factor, it is given by 

Z = {tp\(oo)ip A (x 1 )ip A (x 2 )) = lim x 2hx (jpx(x)tpA(x 1 )ipA(x2)) = (xi - x 2 ) A , (5.5) 

X— ¥oo 

where the exponent is A = h\ 3 — 2h\. The fusion rules (j5.3|) (see also Fig. [5]) indicate that 
h\ 3 = /i2A (k > 1) or h\ 3 = (k > 1). Then using the relation (|4.12p . one arrives at 
A = l/(2(k + 2)) (k > 1) or A = -3/(2(fc + 2)) (k > 1). Correspondingly the partition 
functions are 

Z = (xi -x 2 )~ 1 ^ (k>l), Z 2 = (xi-x 2 ) 1 ^ (k>l). (5.6) 

By inserting them into (|4.17|) and (|4.18|) . the driving processes dxu and dx 2 t characterizing 
the geometric aspects of the SLE traces become: 

a i au , 2a 2 + KAai 2a x + kAo 2 , . 

dxu = Ja\KdB\t H dt, dx 2t = Ja 2 KdB 2t H dt, (5.7) 

Xlt - x 2t " x 2t - X\ t 
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S N 

«A /lA 

/ \ / \ / 

^0 ^0 ^0 

ISLE 2SLE 3SLE 4SLE 

Figure 5: The fusion rules of the bee operator The primary fields with weight h\A f or 
A > k do not exist due to the fusion rules. For instance, for k = 2, the primary fields in the 
shadow area vanish. 



where we normalized the variances by dq a t = a a dt so that d£ at = ^/a^dB a t with two inde- 
pendent standard Brownian motions: E[ciB Qt ] = and E[dB at dBpt\ = 5 a /3dt. These driving 
processes with the SLE (|4.17p describe two curves emerging from two points x\ = x\q and 
X2,o- Setting y s = x\t — X2t and rescaling the time by ds = k(cl\ + 0,2) dt, we reduce the 
processes to the following Bessel process: 

A + 21 k 

dy s = dB s + '—ds (5.8) 

with the effective dimension d e s = 2A + 4//C+L Substitution of the exponent A and k (|4.13|) 
yields 

{I (k = 1) 2 

2 (k+i) n ,v for h \ 3 = °> dgs = 2 + — — (k > 1) for h Xa = h 2 A- (5.9) 

Recalling that the Bessel process is recurrent (resp. not recurrent) if d e s < 2 (resp. d e s > 2), 
we conclude that the driving processes x a t hit each other with probability 1 for h\ 3 = and 
never hit for h\ z = ^A- The hit of the driving processes means the hit of the tips of the SLE 
traces, and hence the case for h\ 3 = describes a single curve (i.e. single arch) whose end 
points are x\ and X2, while the case for h\ 3 = h^h describes two curves converging to the 
point at 00. The above observation agrees with Conjecture 15.11 

5.2 Triple SLEs and arch (crossing) probabilities 

Let us discuss a more non-trivial example: the triple SLEs (m = 3). The geometric properties 
are characterized by the partition function: 

Z = {ipx 4 (x4)ipx s {x 3 )i/)x 2 (x2)ipx 1 (x 1 )}, (5.10) 

where X4 = 00 and by construction the bec operators (1 < j < m) take values in the 
fundamental representation L\. From the fusion rules (see also Fig. [5]), the conformal weight 
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(a) 




(b) 




xi{= 0) x 2 {= x) x 3 (= 1) 



Figure 6: Configurations for m = 3 and n = 1. There are two (c3 t i = (^) — (q) = 2) 
topologically inequivalent configurations. These configurations are described by the partition 
function involving the 4 bcc operators with the same conformal weights h\. The configuration 
in the left (resp. right) panel is characterized by the partition function Zq 1 (resp. Zc 2 )- 



h\ 4 of the bcc operator ip\ 4 (oo) inserted at 00 is h\ A = h 3 \ (k > 3) or h\ 4 = h\ (k > 1). 
When the level takes its value in the range k > 3 and h\ 4 = h 3 \, there exists the factorized 
partition function with the same exponents l/(2(k + 2)): 



[(x 2 -x 1 )(x 3 -xi)(x 3 -x 2 )}^+V (k>3). 



(5.11) 



According to Conjecture 15.11 we expect that this partition function describes three curves 
converging toward 00. 

The case for h\ 4 = h\ is more interesting, because two topologically inequivalent configu- 
rations do exist (see Fig. [6]) . By conformal mapping f(z) = (z—xi)(x 3 —X4 : )/((z—X4 : )(x 3 —xi)), 
the four points Xj are transformed to the points x\ — >• 0, x% — >• x, x 3 — > 1 and X4 — > 00, where 
x = f(x 2 )- Thus the partition function Z can be expressed as 

z = ((x 2 -x 4 )( Xl -x 3 )y 2h ^z(x), z(x) = (ip A {oo)Mi)Mx)Mty)- (5.12) 

The correlation function Z(x) can be calculated by solving the KZ-equation (|4.2ip [22^ 119] . 
Its explicit form reads 

Z(x) = Z Cl (x) + Z C2 (x) 

with 



Z Cl (x) = F[ ] + 



1 



Z C2 (x) 



,(-) 



+ 



1 



,(+) 



(5.13) 
(5.14) 



where F^'s are expressed by the hyper geometric function 2 F% 



,(-) 



-2/i A 



(1-x) 



/l 2 A-2/lA, 



'2-^1 



1 



1 



k 



-:x 



(+) 



(-) 



x 
1 

jfe" 



(+) 



h 2 A-2/iA^ _ a 
- X l-2ftA(l_ a 

-2x h2A ~ 2hA {l 



\/i 2 A-2/iA 2 J7^ 



fc + 2 ' fc + 2 ' k + 2 ' 
1 3 k + A 



; .r 



k + 2 ' k + 2 ' fc + 2 ' 
/c + 3 jfe + 1 Jfe + 1_ 
fe + 2' jfe + 2' A; + 2' ; 
1 3 2 



fc + 2 ' k + 2 ' fc + 2 ' 



; ;r 



(5.15) 
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x 2 {x 




x 2 {x 



Figure 7: Arch configurations in arbitrary domain D obtained by conformal transformation 
of the configurations in Fig. [6) The probability of the occurrence of the configuration (a) 
(resp. (b)) is exactly the same as that of the occurrence of the configuration (a) (resp. (b)) 
in Fig. [6j 



and the coefficients c± are respectively, given by 

r(2/(fc + 2))r(-2/(fc + 2)) 
- r(i/(* + 2))r(-i/(* + 2))' + " 



r 2 (2/(A; + 2)) 



r(3/(k + 2))T(l/(k + 2)Y 
:). For k 



(5.16) 



Note that Zc x {x) and Zq 2 {x) satisfy the relation Zc 1 (x) = Zq 2 (1 — x). 
holds!, and then the partition functions reduce to simple forms 

Z Cl (x) = x- l /\l - xf'\ Z C2 (x) = x l '\l - x)- 1/2 . 



1, c_ 



(5.17) 



For generic k, the behaviors of the partition functions close to the position x = and x 
are described as 



z C x {x) 



-2h A 



(1 



\h 2A -2h A 



X 
X 



Z C2 (x) 



X 



h 2 A-2h A 



(1-x) 



-2/i/ 



X 
X 



(5.18) 



Thus the fusion rules (|5.3|) and Coniecture l5.1 l imply that Zq x describes two curves connecting 
the points [Ox] and [loo] (configuration Ci; see the left panel in Fig[6|) while Zq 2 corresponds 
to the configuration [ooO] and [xl] (configuration C2; see the right panel in Fig[U]). 

Utilizing the partition functions (|5.14p . we can exactly compute the probability P[Ci] 
(resp. P[C2]) of the occurrence of the configuration Ci (resp. C2) with the initial condition 
that the three curves emerging from the point x\ = 0, X2 = x, X3 = 1 by 



P[Ci 



Z Cl (x) + Zq 2 (x) 



P[C 5 



Z C2 {x) 



Z Cl (x) + Zc 2 {x) 



(5.19) 



Since the measure is conformally invariant, the probabilities in arbitrary (simply connected) 
domain D can be easily reproduced by use of Z(x) and the conformal map f^ 1 (z). More 
explicitly, the probability P[Ci] (resp. P[C2]) is exactly the same as that of the occurrence 
of the configuration with two curves joining the boundary points \x\X2{x)] and [X3X4] (resp. 



3 Thus the second terms in (|5.14[1 vanish. In fact, this is a consequence that primary fields with weight h 2 A 
do not exist for k = 1. 
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Figure 8: The crossing probabilities P[Ci] for various level k. 



[xix^\ and [^2(^)^3]) on &D (see Fig. [7]) (note that X2 depends on x, i.e. X2(x) = f^ 1 (x)). In 
this sense, the arch probabilities can be interpreted as crossing probabilities of the interfaces. 
For k = 1, the formula becomes very simple: 



In Fig. [SJ the arch probabilities P[Ci] are depicted for several values of k. As shown in Fig. [51 
the x-dependence of the arch probabilities close to 1/2 in the wide range of x, as increasing 
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